One-Dimensional Impenetrable Anyons in Thermal Equilibrium. I. 
Anyonic Generalization of Lenard's Formula 
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We have obtained an expansion of the reduced density matrices (or, equivalently, correlation 
functions of the fields) of impenetrable one-dimensional anyons in terms of the reduced density 
matrices of fermions using the mapping between anyon and fermion wavefunctions. This is the 
' generalization to anyonic statistics of the result obtained by A. Lenard [l[ for bosons. In the case 

' of impenetrable but otherwise free anyons with statistics parameter k, the anyonic reduced density 

matrices in the grand canonical ensemble is expressed as Fredholm minors of the integral operator 
(1 — 7^t) with complex statistics-dependent coefficient 7 = (1 -f e**'^'')/7r. For k = we recover 
' the bosonic case of Lenard 7 = 2/7r. Due to nonconservation of parity, the anyonic field correlators 

^ , {^1^{x')'^a{x)) are different depending on the sign of x' — x. 

B 

> 

C/3 



I. INTRODUCTION 



^ I During the last few years, several models of one-dimensional anyons [1, H, 0, H, Q have attracted considerable atten- 
^ tion. This attention is motivated, besides general fundamental physics interest of anyons, by several new statistical- 
' , mechanical features of the anyonic systems. One of them is the behavior of the zero-temperature field correlation 
^ functions, which exhibit oscillations with the period dependent on the statistics parameter in the leading term of 
Q their algebraically decaying large-distance asymptotics. These oscillations were obtained within the harmonic fluid 
, approach [?! and conformal field theory 8] . 

' Much of this effort was focused on the anyonic model proposed in [1] and clarified in Q . This model can be 
understood as an anyonic extension of the Bose gas with 5-function interaction [§] which was solved by Lieb and 
^ . Liniger. The study of this model was initiated by Batchelor et al. [H!, [H, and continued in 0, [l,[i3l- In the hmit 
f^^ — of infinite interaction strength, the anyons become impenetrable. The present work is the first paper in a series that 
0^ intends to provide a complete treatment of the correlation functions of such impenetrable free anyons in a manner 
similar to the impenetrable Bose gas [Tsj . The first step taken in this work is the anyonic generalization of Lenard's 
, formula [l| which gives the reduced density matrices of hard-core bosons in terms of the fermionic reduced density 
^— I matrices. With this result, the reduced density matrices of free impenetrable bosons can be expressed in terms of 
Fredholm minors of the integral operator Ot with the kernel given by the Fourier transform of the Fermi distribution 
QP function. 

In this work, wc use the Anyon-Fermi Q and Anyon-Bose ^3] mapping to generalize the treatment of Lenard 
to the case of anyons. Our main result is Theorem IIV.II which in the case of impenetrable free anyons leads to a 
representation of the correlation functions of anyonic fields as Fredholm minors of the operator Ot. This representation 
of the correlation functions can be summarized as follows. The Fermi distribution function ■d{k,T, h) at temperature 



5^ , T and chemical potential h defines the integral operator 6t with the kernel 

1 f°° 1 

ST{x^v) = -j dke'^'^--y^d{k,T,h), d{k,T,h)= (1) 
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which acts on an arbitrary function f{x) as 

(eTf){x)^ jeT{x~y)f{v)dv, (2) 

where / is an interval or a finite union of intervals on the real axis. The resolvent kernel qt{x, y) associated with the 
kernel 9T{x,y) is defined to satisfy the equation 

Qrix, y) - 1 I 9t{x - z)qt{z, y) dz ^ Orix - y) . (3) 



Then, in the thermodynamic limit, the temperature-dependent correlator of the anyonic field operators is given by 
the following relations: 

• If a;' > X, 

{^\{x')-^A{x))T,h = -Qt{x', x) det (l - 7^t) , (4) 

IT \ / 7=(l+e+""')/7r 

where det ~ "f^rj is the Fredholm determinant of the integral operator 9t (see Appendix [X]) which acts on 
the interval /+ = [a;,a;']. 

• If x' < X, 

{^\{x')^Aix))T,h = -gTix,x')det(l-jeT) , . , (5) 

and now the integral operator acts on the interval /_ — [x',x]. 

Using Eq. (|A11[) we can also express these formulae in terms of the first Fredholm minor of the integral operator 6t- 
The methods used in our work allows for direct extension of these results to the 2n-point correlators with certain 
ordering of the arguments: 

• If xi < x[ < ■ ■ ■ < Xn < x'„, 

{^^AO ■ ■ ■ ^l{x[)^Aix^) ■ ■ ■ ^a{x^))t,h = : ' """^ gT f g' • • • ' ) det (l - 7^t) ^ ^. ^ , 

(6) 

where now the integral operator acts in the region = [xi, x'l] U • • • U [a;„, x^], the statistics factor C is defined 
by Eqs. ([?T|) . (|^^ . and ([^51) . and other notations are given in Appendix \K\ 

• li x'l < Xi < ■ ■ ■ < x',^ < Xn, 

{^lix'J ■ ■ ■ ^i{x[)^A{x^) ■ ■ ■ ^AiXn))T.H ^ ' ' "^"^ gT f g ' ' ' ' ' g ) det (l - ^er) ^ . ^ , 

TT" \Xi,- tX^J \ J 'y=(l+e-""<)/7r 

(7) 

where the integral operator acts in the region J_ = \x'y^ x\\ U • • • U [x^, a;„]. 
Again, Eq. (jAlip can be used to express the 2n-point correlators as the nth Fredholm minors of the integral operator 

The paper is organized as follows. In Sec. |TT1 we introduce the reduced density matrices of anyons and their relation 
to the correlation functions of the anyonic fields. In Sec. IIIIl we describe the mappings between the wavefunctions 
of the one-dimensional particles with Fermi, Bose, and Anyonic statistics. In Sec. IIVI we prove the main theorem 
which expresses the reduced density of anyons in terms of the reduced density matrices of bosons and fermions. The 
consequences of this theorem in the case of free impenetrable particles is analyzed in Sec. |Vl where we obtain the 
reduced density matrices of anyons as Fredholm minors of an integral operator. The basic information on Fredholm 
determinants, details of calculations with anyonic fields, and a sketch of the derivation of the reduced density matrices 
for free fermions in the grand canonical ensemble are presented in three appendices. 



3 



II. FROM CORRELATION FUNCTIONS TO REDUCED DENSITY MATRICES 

The one-dimensional anyons considered in this work are characterized by anyonic fields ^'^(a;), ^'^^(a;) which obey 
the following commutation relations 

^a{xi)^\{x2) = e-'^'''^^'-^^^^\{x2)^A{xi) + S{x, - X2) , (8) 

^a{xi)^a{x2) = e'^-<=''-^-'>^A{x2)'i'A{xi) . (10) 

Here k is the statistics parameter, and e{x) = x/\x\, e(0) = 0. The commutation relations become bosonic for k = 
and fermionic for k = 1. For an arbitrary Hamiltonian of the anyons confined to the interval V = [—L/2, L/2\, the 
TV-particle eigenstates are defined as 

i^Niim = / • • • / dz^ x%{zi, z^\{\})^>\{zm) ■ ■ ■ ^\{zim , (11) 

ViV! Jv JV 

(1'jv({A})| - -L= ( dzi-- - [ dzN {0\^a{zi) ■ ■ ■ ^a{zn)Xn{zi, • • • , zn\{X}) , (12) 



/TV! 



V 



where Xn the (norm one) quantum-mechanical wavefunctions of N anyons, and {A} is a set of quantum numbers 
specifying the state. Note the order of the field operators in these relations, which is dictated by the boundary 
conditions on the wavefunctions in the periodic or quasi-periodic case @. The wavefunctions x have the anyonic 
symmetry 

XAr(2i, • • • ,Zi,z,+i, ■■■ ,zn)^ e"'*'(''*~^'+^^X?^r(2i, • • • ,Zi+i,z,, ■■■ ,zn), (13) 

that reflects the field commutation relations. We are interested in computing the finite-temperature correlation 
functions of anyonic fields. The simplest example of these correlators is 

{^iix')^Aix)). (14) 

In the grand canonical ensemble characterized by temperature T and chemical potential h, the field correlation 
function is given by following relation 

°° p-E{{X})/T 

{^\{x')^a{x))t,, = ^^e"^/^ {^^{{Xml{x')^A{x)\^Ni{X})) , (15) 

^=1 {A} { y y ) 

where E{{X}) is the energy of the eigenstate with quantum numbers {A}, and Z{h,V,T) is the grand-canonical 
partition function 

oo 

Z{h, V,T)^Y.J2 e'^^/^e-^"^"/^ . (16) 

N=0 {A} 

As shown in the Appendix [B] the correlator at a fixed number of particles N is given by an overlap integral of the 
corresponding wavefunction 

{^N{{X})\^\{x')^Aix)\'^Ni{X})) =N [ dzi--- [ dzN^lXNi^U--- ,ZN-l,x'\{X})x%{zi,--- ,ZN-1,X\{X}), (17) 

Jv Jv 

so that the correlation function (fT5|) can be written as 

{¥^ix')^Aix))TM = ^^e^^/^ " N dz,--- dz^., 

N=l {A} [ y y ) Jv Jv 

xxn{zi,--- ,ZN^i,x'\{X})x%{zi,--- ,ZN-i,x\{X}). (18) 
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We will also be interested in a class of 2n-point field correlation functions at finite temperature: 

(vf^K) • • • ^\{x',)^a{x,) ■ ■ ■ ^A{Xn))T,k ■ (19) 

These correlators can be expressed similarly to Eq. ([T51l 

°° p-E{{\})/T ATI /■ r 

(vl/t^(x;)...vI^^(x;)M>^(xi)...^^(x„))T..= E^e^'^/^ ^^^^^^^^ y^dzi... y^dz^_„ 



|{A}). (20) 

As one can see from ([T51) and ([201) , the correlation functions are obtained as a combination of the wavefunctions and 
ensemble probabilities. In this context, it is useful, similarly to the case of fermionic or bosonic particles, to introduce 
the reduced density matrices of anyons: 

Definition 1. For a statistical ensemble characterized by the probabilities pf\j, the anyonic n-particle reduced density 
matrix is defined as 

(Xl,--- , Xn\pl\x[, ■ ■ ■ ,x'^)= ^ ^Vfx} y / dzi--- f dzN-n 

N=n {A} ^ •'^ •'^ 

X Xn{zi,--- ,ZAr_„,a;i,- • • ,x'J{\})x%{zi, ■ ■ ■ , ZAr_„, Xi, • • • ,a:„|{A}), (21) 
where the wavefunctions x% {a} '^'"^ normalized to one. 

If the probabilities coincide with those in the grand canonical ensemble, p^-^y — e^'^^'^e^^^-^^^^^'^ /Z{h, V,T), 
the one-particle reduced density matrix is just the 2-point correlator (jlSp . and the n-particle reduced density matrix 
is the particular 2n-point correlator (j20p . These relations are exactly the same as in the case of bosonic and fermionic 
statistics. A particular "anyonic" feature of the Definition [Tl is the fact that we integrate over the first N — n 
arguments of the wavefunctions. In the case of bosonic and fermionic reduced density matrices, integration over any 
subset of the N — n out of N arguments produces the same result due to the parity of the wavefunctions. This is not 
the case for the reduced density matrices of anyons due to the anyonic symmetry p3p which in general, e.g., in the 
periodic or quasi-periodic situation ( "anyons on a ring" ) makes different arguments of the wavefunctions inequivalent 
- see discussion in the next Section. 

Under a certain set of conditions, which is also made precise in the next Section, there is a correspondence between 
the anyonic and the fermionic or bosonic wavefunctions. This correspondence will be used later to express the anyonic 
reduced density matrices as expansions in terms of fermionic or bosonic ones. 



III. ANYON-FERMI AND ANYON-BOSE MAPPING 



To establish the correspondence between the wavefunctions of anyons and fermions or bosons we define the two 
functions which essentially incorporate statistical properties of the wavefunctions of different statistics in one dimen- 
sion: 

A^izi,--- ,z^)=e-''5:,<.^(^.-^.)/2 (22) 

and 

B{zi,--- ,ZN)^Y[e{zj-Zk), (23) 
j>k 

where the notations are the same as in Eq. ([5]). 

The mapping between anyons and fermions or bosons is analogous to the Bose-Fermi mapping discovered in p^ . 
where it was noticed that any wavefunction of N fermions has a bosonic counterpart given by 



x''{zi,--- ,zn) =B{zi,--- ,ZN)xH^ir-- ,zn) 



(24) 
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This correspondence is valid under very general conditions, with no restrictions on the external or particle-particle 
interaction potential, except for the requirement of the hard-core condition which should make the bosons impenetra- 
ble. For particles confined to a box with "hard wall" boundary conditions (BC), bosonic and fermionic wavefunctions 
satisfy the same BC. In this case, if is an eigenfunction of the Hamiltonian, then x** is also an eigenfunction with 
the same eigenvalue. However, in the case of a ring of length L with periodic BC for bosons, the BC for the fermions 
are in general different and given by 

X^(0,-.. ,ZAr)-(-l)^-V(ir- - ,^iv). (25) 

For even TV, when Eq. (P5|) means that the BC for fermions and bosons are different by a phase shift tt, the relation 
between the eigenenergies of the fermionic and bosonic systems is less direct. Since for non-coincident coordinates 
= 1, the Bose-Fermi mapping (|24p is symmetric and remains true if the superscripts b and / are interchanged. 



A. Anyon-Fermi Mapping 

It is straightforward to see that similarly to the Bose-Fermi mapping ([M)) . the wavefunction with anyonic symmetry 
(fT3|l can be obtained by multiplication of a fermionic wavefunction with the statistics factors ([22]) and ((23|) 0, 0] : 

X^izi,--- ,zn) ^A^,{zi,--- ,zn)B{zi,--- ,zn)xHzi,--- ,zn), (26) 

Besides the anyonic symmetry (IT3l) . the wavefunction x° ([26| satisfies the condition 

X'^(zi,--- ,Z7v)|.,=., -0 foraU{i,j}e{l,... ,iV}. (27) 

This means that the correspondence (j26p is valid as long as potential energy contains a hard-core part which ensures 
that the anyons are impenetrable and condition (|27|) is indeed satisfied. Other properties of the Anyon-Fermi mapping 
(|26p are similar to those of the Bose-Fermi mapping. It is valid for an arbitrary form of the potential energy in the 
particle Hamiltonian. When the particles are confined to an interval with "hard wall" boundary conditions, the 
fermionic and anyonic systems both satisfy the same BC of the wavefunctions vanishing at the ends of the interval. 
In this case, if is an eigenfunction of the particle Hamiltonian, then x° is also an eigenfunction of this Hamiltonian 
with the same eigenvalue. This follows from the fact that the statistics factors (|22|) and ((23|) are constant everywhere 
except for the points of coincident coordinates, where the wavefunctions vanish. 

When particles are confined to a ring with periodic or quasi-periodic BC, the properties of the Anyon-Fermi mapping 
are more complicated. In this case, the anyonic wavefunction will have different boundary conditions for each of its 
coordinate (see Q and Appendix A of Q), the difference being given by an extra phase shift that depends on the 
statistics parameter k. Specifically, if the fermion wave function obeys some generic quasi-periodic BC (the same in 
all coordinates) which can be written as 

X^(0, ■■■,zn)^ {-lf-'e-^*xHL, ■■■,zn), (28) 
then the anyonic wavefunction obeys the following BC in its different arguments 

X'^(0,Z2,-- - ,ziv) = e-'^ x"(i,^2,--- ,ziv), 
X'^(zi,0,--- ,zw) = e'(2— «x"(^i,^,--- ,^iv), 

: _ (29) 

X"(zi,z2,--- ,0) = e^(2(^-i)-''-^)x'^(^i,^2--- ,i), 

where cj) = (j)-\-TTK{N~l). As for the Bose-Fermi mapping ([M)) with even N, the anyonic and fermionic eigenenergies are 
not related directly in the situation of a ring with quasiperiodic BC. In physics terms, this difference between anyons 
and fermions corresponds to the statistical magnetic flux ttk^N ~ 1) through the ring produced by N one-dimensional 
anyons of statistics k. 

The Anyon-Fermi mapping (I26p is not symmetric. The inverse relation can be written as 



X {zi, ■ ■ ■ ,zn) = A^k{zi, • • • , zm)B{zi, ■ ■ ■ ,ZAr)x°(2;i, ■ • • ,zm) 



(30) 
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B. Anyon-Bose Mappping 

The Anyon-Bose mapping was historically the first mapping of this kind introduced for one-dimensional anyons in 

X"(zi, • • • ,zn)^ A^{zi, ■ ■ ■ ,zn)x''(zi, ■■■ ,zjv) . (31) 

The wavefunction x° in (|3ip has the correct anyonic symmetry (|13p . and in contrast to Anyon- Fermi mapping (j26p . 
need not vanish when a pair of coordinates coincide. It should be noted, however, that without this hard-core 
condition, the discontinuity of the statistics factor (|22p at coinciding coordinates translates into discontinuity of 
the wavefunctions l|3ip . In this case one needs an additional condition regularizing the wavefunctions. Also, without 
the hard-core condition, the statistics factor A changes substantially the behavior of the wavefunctions at the points 
of coincident coordinates (for instance, if the particle-particle interaction is (5-functional, statistics renormalizes the 
interaction strength Q) and the energy eigenvalues of the bosonic and anyonic problems are different regardless of 
the boundary conditions. 

Other properties of the Anyon-Bose mapping (j3ip are very similar to those of the Anyon- Fermi mapping. It is valid 
for arbitrary potential energy. For particles in a box with "hard wall" BC, both wavefunctions (I3ip satisfy the same 
condition = and = at the boundary. For particles on a ring with generic quasi-periodic BC for bosons that 
can be written as 

X^O,--- ,z^) = e-^V(i,-- - ,^iv) (32) 
(and have the same form for all other arguments of x''), the anyonic wavefunction obeys the following BC 

X'^(0,Z2,-- - ,;2iv) = e-'^ x"(i,^2,--- ,ZAr), 
X"(zi,0,--- ,zjv) = e'^'^--^^x''izi,L,--- ,zn), 

'■■ _ (33) 

X"(Z1,Z2,--- ,0) = e^(2(7V-l)..-0)^a(^^^^^... 

where (j) = cj) + ■nK{N — 1). The Anyon-Bose mapping is also not symmetric. The inverse of ([3T|) is 

X^{zi, ■■■ ,zn) = A_k(zi, • • • ,zj^)x°'{zi, ■■■ ,zn)- (34) 



IV. MAIN THEOREM 



In this section, we consider an arbitrary statistical ensemble in which the states Xn {a} occur with probabilities 
pfxy The anyons are assumed to be confined to an interval V = [—L/2, L/2], and wavefunctions are normalized 
to 1: ||xAr{A}ll — 1- ^'^^ §0^1 estabhsh a relation between the reduced density matrices p° of anyons (^J) 

and similarly defined reduced density matrices of bosons and fermions p^. The fermionic and bosonic states 
that correspond to the anyonic states x^ in the Anyon- Fermi (j26p or Anyon-Bose (|3ip mappings have similarly 

normalized wavefunctions, and we assume that they have the same probabilities pfx}- This assumption is natural 
under the conditions (discussed in the previous Section) for which the energies of the states of different statistics are 
the same, as they are, for instance, when the wavefunctions satisfy the "hard wall" boundary conditions and hard-core 
condition on the particle-particle interaction. The relation between the reduced density matrices is established by the 
following theorem: 

Theorem IV. 1. Let xi, • • • , a;„, x']^, • • • , be 2n coordinates in the interval V , and 0± are the parts of the space 
of these coordinates in which they are ordered, respectively, as xi < x'l < ■ ■ ■ < Xn < x'^ and x[ < xi < ■ ■ ■ < 
x'^ < Xn- For the 0+ ordering, one can define the subset of V: /+ — [xi,x'i\ U • • • U [a;„,a;^] C V, and the subset 
/_ — [x'l , xi] U • • • U [a;^, Xn] C V for ordering as in 0_ . // the conditions of validity of the Anyon-Fermi I126\) or the 
Anyon-Bose mapping i31\) are fulfilled, the reduced density matrices of anyons can be expressed then in terms of the 
reduced density matrices of fermions as 

(xi, • • • , • • • , x'j^)± = A^k{xi, • • • , x'j^)B{xi, • • • , x'j^)Ak{xi, • • • , Xn)B(xi, • • • , a;„) 
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X V(-l)-' / dzi--- dzj {xi,- ■ ■ ,Xn,zi,- ■ ■ ,Zj\pi_^_Ax[,- ■ ■ ■ ■ ,Zj) (35) 

or bosons as 

{Xi,--- ,Xn\pl\x[,--- ,x'J± =A_i^{x[,--- ,x'JA^{xi,- ■ ■ ,Xri) 

xV(-l)-' / dzi--- dzj (xi,--- ,Xn,zi,--- ,Zj\p'',^^j\x[,--- ,x'^,zi,--- ,Zj). (36) 

,-n 3- Jl± Jl± 



j=0 



The subscript ± in these expressions specifies whether Xi, • • • , Xn, x[, ■ ■ ■ , x'j are ordered as in or 0_. 

Proof. The proof follows that of Lenard [IJ, generalizing it to the anyonic statistics. As a first step, we need a 
preliminary result. 

Lemma IV. 1. For any symmetric function f{zi, • • • , z„) and a constant a, 



I dzi ■ ■ 


■ 1 dzj / dzj^i ■ ■ 


■ 1 dZnfizi, ■ ■ 


I ) 7 


'i± 


Jl± JV 


JV 





j=0 

where C" = {n-jy.j\ '^"■'^ number of variables zi, - ■ ■ , Zn contained in I±. 

Proof. The L.H.S. of ([37]) can be written explicitly as 



'v\i± 



/ dzi ■ ■ 


■ / dz,n 1 


'i± 


J I± JV 



(38) 



and combined with an obvious relation /y\/^ dzi = Jy dzi — Jj^ dzi , {i = m + 1, ■ ■ ■ ,n) can be further transformed 
into 



g = V C;^,a" "V <^r™(-l)' / dzi--- f dz,n+k I dz^+k+1 ■ I dznfizi, • • • , z„) . (39) 

m=0 k=0 •''± ''I± ''V 

Collecting the terms in this expression with the same j = m + fc, we obtain the desired result 



3=0 



/ dzi- ■ 


■ / dzj / dzj+i • • 


■ / dZnf{zi, ■ ■ 


I ^n) 




J/± JV 


JV 





□ 



Now we can prove the Theorem lIV.il starting with (|35p . Using the Anyon- Fermi mapping (|26p we have 

{Xi,--- ,Xn\Pn\x'ir ■ ■ ,x'^)± = ^ ^PW (]\jJ )\ / ^.Zi ■ ■ ■ / dzN ^^C {xi, ■ ■ ■ , Xn, x[, . . . , x'„)± 

N=n {A} ^ ^ ^ 

X Xjv(^l'-" ,2Ar-n,a;i,--- , x'J{X})Xn{zi, ■ ■ ■ ,ZN-n,Xi,--- ,X„\{X}), (40) 



where 



C'{xi,'-' ,Xn,Xi,-'' ,x^)± — A—f^[xi,-'' , Xjj)i3(x]^ , • • ■ , 2;„) (xi , • ■ ■ , a;n)i3(a;i , • • • ,2;^) 

n N~n 
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One can see directly that 

j = l 

This means that 

C{xi,- ■■ ,Xn,x[,--- ,x'^)± ^ • • • ,x'n)B{x[,- ■ ■ ,2:^)71^(2:1, • • • ,Xn)B{xi, ■ ■ ■ , a;„)(-e^"'')'^ , (43) 

where a'{I±) is the number of variables zi, • • ■ , zj^^n in I±- Applying now Lemma llV. II with a — — e^*'^'*, we obtain 
for the anyonic reduced density matrices 

°° TV! 

{xi, ■ ■ ■ ,Xn\PnWi^ ■ ■ • ,x'j^)± = A^f^{x'i, ■ ■ ■ ,x'j^)B{x'i, ■ ■ ■ , a;^)v4„ (^i , • • • ,2;„)_B(2:i, • • • ,Xn) ^ Pu> / lyr — 

N-n r n n r 



X }^C;-"(-l)^(l + e±"T / dz,--- dz, / / x:v,{A}X^,{A} • (^4) 

Interchanging the order of summations, one can notice that the sum over A'' and {A} is precisely p^j^j- Therefore 
finally 

{X\, ■ ■ ■ ,Xn\Pn\Xi,''- ,Xn)± = A-k{Xi, ■ ■ ■ , X^)B[Xi, ■ ■ ■ , X^)AKixi, ■ ■ ■ , Xn)B(xi, ■ ■ ■ , Xn) 

X V'(-l)^ r. / dzi--- dzj {xi,- ■ ■ ,Xn,zi,- ■ ■ ,Zj\pI_^_Ax[,- ■ ■ ,x'^,zi,- ■ ■ ,Zj) . (45) 

The proof of (|36p is similar. In this case, we use the Anyon-Bose mapping (pij) . and the C± function is 

C{xi^--' ^XmX-^^-- - ,2:^)^ = A — ^ix-^^-- - , 27^)yl^ (xi , ■ ■ • ^Xri) 

n N — n 

■J~J ^-™Ke(zi-x'j)/2^iTrKe{zi~Xj)/2 ^ ^^g-j 

3 = 1 1=1 

This means that we can use Lemma llV. 1 1 with a = e**'^". Interchanging the order of summation and identifying the 
bosonic reduced density matrices p^^^j we obtain (|36p . □ 

The results of Theorem IIV.II do not depend on the statistical ensemble used in the computation of the reduced 
density matrices as long as the particles are subject to the hard-wall boundary conditions making the state energies 
independent of the statistics. They also do not depend on the form of the interparticle potential beyond the need for 
the hard-core part which ensures that the wavefunctions satisfy the hard-core condition. If both of these conditions 
are satisfied, we can see from (|35|) and ((36)) that there is also no explicit dependence on the length L of the confining 
box V , and the results remain valid in the thermodynamic limit L ^ 00. 



V. IMPENETRABLE FREE CASE 



The Anyon-Fermi relation derived above for the reduced density matrices is particularly useful in the situation 
when the radius of the hard-core interaction is vanishingly small, and no other interactions are present. In this case, 
the fermionic problem is identical to free fermions, since the hard-core potential of zero radius effectively vanished 
due to antisymmetry of the wavefunctions. The reduced density matrices p(^ coincide then with those of free fermions 
[H (see Appendix [C|: 

(xi, • ■ • , xMWi, ■■■^0-^(^t(1}'["'17), (47) 
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where OT{x,y) /it is the Fourier transform of the Fermi distribution function: 



dk- 



Jk{x-y) 



At T = we have 



1 + g(fc2-ft)/T 

sin q{x — y) 
x-y 



(48) 



(49) 



where q — \fh is the Fermi momentum. 
Applying Theorem IIV.II and (jC8|l we have 



(ii,--- ,x„|p^|a;'i,-- - ,0± = C^(-1)- 



(1 + e±"«)^' 



i± 



where 



C(a;]^ , • • • , Xn) = A_^(a;-^ , • • • , a;„)_B(a;j^ , • • • , x^^) 71^(0:1 , • • • , a;„)-B(a;i , • • • , 



(50) 



(51) 



and the subscript ± specifies particular ordering of a:i, • • • , x„, a;^ , • • • ,x'^ as in Theorem IIV.II This result can be 
rewritten in terms of Fredholm minors using (jA10|) 



(xi,--- ,Xr,\p'i^x\,--- ,x'^)± = —C{x[,--- ,a;„)det 1-76*^ 



-^1 7 ' ' ' J Xn 



7 = (l+e±''"=)/7r 



where the integral operator 9j, with kernel 9T{x,y) is defined by its action on an arbitrary function /: 



-{x,y)f{y)dy 



i± 



Finally, introducing the resolvent kernel g^{x,y) associated with the the kernel Orix^y), which satisfies 

(1 + e±"'^) 



QT{x,y) 



and making use of (|A11[) . (|^D|) and (US]), we obtain 



-(x - z)gj.{z, y)dz = Orix - y) , 



{^\{x'^) ■ ■ ■ ■^\{x',)^a{xi) ■ ■ ■ ^A{Xn))T,h.. 



C(Xi, • ■ ■ ,Xn) ± ( Xi, - ■■ ,Xn 



Qt 



det 



7 = (l+e±''"=)/7r 



In the particular case of the simplest two-point correlator, this expression reduces to 

{^\{x')'i'A{x))T,h.,± = ^QTi^'^x) det (1 - 7^±) 



7=(l+e±"'=)/7r 



(52) 



(53) 



(54) 



(55) 



(56) 



and gives the correlator of two anyonic fields in terms of the Fredholm determinant of the integral operator 9t and 
its resolvent kernel. 



VI. CONCLUSIONS 



In this work, we have studied the field correlation functions (reduced density matrices) for free impenetrable anyons 
in one dimension and obtained the representation for these functions in terms of the Fredholm minors of the integral 
operator with the kernel given by (j48p . This representation of the correlation functions generalizes similar results 
[1] for one-dimensional impenetrable bosons, and can be used to study the asymptotic behavior of the correlation 
functions beyond the approximation of conformal invariance. In the case of bosons, there is also an alternative 
representation of the zero-temperature one-particle reduced density matrix as the determinant of a Toeplitz matrix, 
development of which was motivated by the study of momentum distribution of bosons in the ground state [l^, [H, [13 ■ 
This representation was extended recently to anyons by Santachiara et al. U in their study of the entTngtat 
entropy of impenetrable free anyons. As in the case of bosons, explicit demonstration of the equivalence of the two 
representations for anyonic correlation functions is an open problem. 
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APPENDIX A: FREDHOLM DETERMINANTS 

In this Appendix, we give a brief summary of results of Fredholm theory of integral equations. For more details, 
see, e.g., [1^. Consider the Fredholm equation of the second kind 

f{x)--i ( K{x,y)f{y)dy = g{x), (Al) 

J a 

where the kernel K{x, y) is a symmetric, bounded and continuous function. 

Defining operations with kernel K{x, y) similarly to the usual matrix operations: 

K^{x,y)= ( K{x,z)K^-\z,y) dz, with K\x,y) ^ K{x,y) , (A2) 

J a 

and 



TtK = / K{x,x)dx, TtK'^ = / / K{x,y)K{y,x) dxdy , and so on, 

J a J a J a 



(A3) 



we have the formulae that are useful for calculation of the Fredholm determinant of the integral operator 1 — ^K: 

(l-7ii:)"i = 1 + 7^1 +72A'2 + --- , (A4) 

and 

lndet(l-7K) = -^^TrX". (A5) 



Indeed, writing (|A5p as 

det(l - 7^:) = IT cxp{-^Tr if"} , (A6) 

n=l 

and collecting terms of the same order in 7 one can see that the determinant can be written conveniently as 

det (1 - ,k) . E(-l)"^ X ■ ■ la 

where 

KJlY--/lA^ dei [K{x,,y,)]. (A8) 

\ yl: : tin J l<j,k<n 

The resolvent kernel R(x, y) associated with the kernel K{x, y) is defined as i? = (1 — 'yK)~^K, i.e., 

R{x,y)--f [ K{x,z)R{z,y)dz = K{x,y). (A9) 



If one introduces the determinants i?„ of kernels R similarly to (lASp . an important relation can be proven to exist 
between i?„ and the r-th Fredholm minor defined as a natural generalization of Eq. \A7\ : 

(AlO) 



det ^1 - -fK 
The relation is Ui 



!:0-E(-i)"4/ dx,---( dx^K^ (y^^ • 



Vi, ■■■ , Vr J ^ n\ J„ J„ ' V yi' ■ ■ • ' Vr^ • • • , Xr. 

n— 
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APPENDIX B: ANYONIC CORRELATORS 



In this Appendix, we prove Eq. (|17p following the approach used in ^ for calculation of the anyonic matrix elements. 
We start with the simple case of the correlator 

{^2\'fi{x')^A{x)\^2), (Bl) 

where we have omitted the quantum numbers {A} unimportant for the present computation. From (jlip and (|12p . we 
have 

(^2\^^^ix')^Aix)\^2) = \j dz'dy' xr(yi,2/2)x5(^i,^2)(0|*A(yi)*A(2/2)*Ii(a;')*A(a;)vI'Ii(z2)vI/Ii(zi)|0). (B2) 
Defining for the moment 

A = mA{yi)^A{y2)^\{x')^A{x)^\{z2)^\{zim , (B3) 
and using the commutation relation ([8]), ^'^(2^)10) — 0, and (0|0) = 1 we obtain 

A = {0\^A{yi)'i'A{y2)^i{x') [^\{z2}^A{x)e-^'^^<^''-^ + S{x - Z2)] ^\{zi)\0) , 

= (0|vI/^(yi)^'A(y2)*!i(^')*!i(^2)*A(a:)*^(zi)|0)e-"''^(^-^^) + ^ A{yi)^ A{y2)'^l{x')^\{z^mSix - Z2) , 
= {Q\^A{yi)-^A{y2)^\{x')^\{z2m 5(x - zi)e-""^(^-^^) + (0|f A(2/i)4'A(2/2)*i(a:')*i(^i)|0) 5{x ~ Z2) . (B4) 

V ' ^ v ' 

(a) (b) 

Performing similar transformations we find that 

a = <S(yi-a;')%2-^2)e-^""^(^^""')+<5(y2-a;')%i-22), 

b = <5(yi - x')8{y2 - z,)e-'^^<y--^'^ + 5{y^ - Zi)5{y2 - x') . (B5) 
Substitution of (|B4| and (|B5]) into (lB2| gives 

{^2\^\{x')^a{x)\^2) dz^[xT{x ,z,)x'i{x,z,)e-^^^^<^^--'^+<^-^^'^^ + xTizu^^^^^^ 

XT{^': zi)x^(zi, a;)e-'--(^^--') + xT{zux)x''2{zux)] . (B6) 

Anyonic property of the wavefunctions (|13p together with its complex conjugate 

X*n{z^: • • • , ■■■ ,zn)^ e-""^("'-"'+i)xw (^1. • • • . z,+i,z,, • • • , z^) , (B7) 

means that (jB6l) reduces to a simple form 

{■^2\^\{x')^a{x)\^2) =2 fdzi xTizi,x')x^^{zi,x) . (B8) 



The generalization to the N-particle eigenstate is straightforward and gives 

{^n\-^\{x')-^a{x)\^n) f dz""-' XNizi, • • • , ZN^,,x')x%{zi, • • • , zjv-i, a;) . (B9) 



APPENDIX C: REDUCED DENSITY MATRICES OF FREE FERMIONS 

The reduced density matrices for free ID fermions were calculated in the original paper of Lenard [l|. To make our 
discussion self-contained, we provide here a sketch of the proof and the main results in the notations that in general 
allow for an external potential U{z) acting on the particles. 

We assume that the fermions are confined to the domain V = [~L/2, L/2], and have a complete set {u\{z)} of 
normalized single-particle wavefunctions with energies e\ in the potential U{z). For instance, for U{z) = 0, and the 
"hard wall" boundary conditions at the boundaries of the domain V, 



, sin(Az), A = , 

Mz)={ \ \ , (CI) 

'fcos(Az), A=f,^,.--, 
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and, with appropriate conventions, = A^. The iV-body wavefunction of a stationary state is given by the Slater 
determinant 



N 



f 

Xn 



(zi,--- ,zn\{X}) = n "A, (^ttw), 



(C2) 



where the set {A} consists of non-coincident single-particle states A^, and the energy eigenvalue is E{{\}) — X^fli 
In the grand canonical ensemble, the Gibbs measure is 



,-E({\})/T 



with Z(h,L,T) = XUI' 



,hN/T~E{{X})/T 



(C3) 



N=Q {A} 



where h is the chemical potential. Using the fact that, with an extra factor 1/A^! included to compensate for over- 
counting, summation over {A} can be replaced with summation over independent individual A^'s, one obtains the 
following fundamental formula 



hN/T-E({\})/T*f 



{A} 



Here we have used (IA8I) and 



X*J{zu--- ,z^\{\})xU^[,--- ,z'^\{\}) = ^ E (-irrin^.,4w) (C4) 

(C5) 



1 P / zi, • • • ,zjv 



(C6) 



Equation (|C4[) and proper normalization of the wavefunctions u\ show that the fcrmionic statistical sum (jC3p can be 
expressed as the determinant (IA7p of the integral operator with kernel (|C6P : 



Z{h, L,T) = y2 ^ j dzi - ■■ j dzN Fn 



Zl, • • • ,ZN 
Zl, • • • ,ZN 



det(l-hF). 



Similarly, using the definition (jAlOp of Fredholm minor of the same operator we see that 



-E({\})/T_ 



Nl 



N= 



{A} 



(N-ny. Jy 



dzi 



dZN-n 



(C7) 



XNizi,--- ,ZN-n,xi,--- ,Xn\{\})xNizi, - ■ ■ ,ZN-n,x[,--- ,a;^|{A}) = dct ( 1 -|- F 









; Xn \ 






Xi, ■ ■ 





SO that the reduced density matrix of the fermions can be expressed as: 

,<) =det (l + F 



{xi, ■ ■ • , a;n|p;^|a;]^, ■ 



^1 , ' ' ' , Xn 
■^1 1 ' ' ■ I ■^n 



/det(l-HF). 



(C8) 



(C9) 



The relation (|A11[) for Fredholm minors means that this result can be expressed simply in terms of the resolvent 
kernel 6Tix,y)/TT associated with kernel F{x,y) (|C6p (factors of tt are chosen so that the notations are the same as 
in the main text): 



{Xl, • • ■ , Xn\Pn\x\, ■ ■ ■ 



I I 



In the thermodynamic limit with no external potential, U{z) = 0, the resolvent kernel 6t is given by 



(CIO) 



(CU) 
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